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In multiquark correlator analyses with 1/N C classifications, it is possible to separate 
the scattering background and to justify the factorization of condensates, which 
allows us to achieve an isolated peak saturation in the QCD sum rules for multiquark 
currents. Then we can extract leading N c properties of the ground state. An 
application to the a meson is demonstrated. 



1 QCD sum rules for multiquark correlators: Difficulties 

The multiquark states have been intriguing subjects of researches to study the infrared 
aspects of QCD beyond those manifested in the usual hadron spectra. One of the con- 
ventional approaches from QCD is to analyze correlators of multiquark currents, which 
is related to the hadronic spectral functions through the dispersion relation. Compared 
to analyses for usual meson and baryon correlators, the multiquark cases are expected 
to have the following difficulties: (i) The complications associated with the many-body 
system, which make the calculation of correlators much involved, (ii) The decays into 
usual hadrons, whose backgrounds contaminate the signal from multiquark states. We 
will show that the use of 1/N C expansion [TJ [5] greatly reduces both difficulties, and, 
more importantly, enables to make logically clear statements [3]. 

We perform the correlator analyses using the Borel transformed QCD sum rules 
(QSR) for the hadronic correlators J dxe lqx (T 'J(x)j1 '(0)) 0]: 
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where LHS is calculated by the operator product expansion (OPE) and RHS is the 
integral of the hadronic spectrum. Here M characterizes the scale of momentum flow 
between interpolating fields, and [i is factorization scale separating high energy calculable 
part and low energy condenstate part. 

Usually we are interested in the ground state properties of RHS in eq. (JTJ) , and they 
can be investigated by taking small M 2 due to the exponential factor e~ s l M ~ . But we can 
not take arbitrarily small M 2 value since the OPE is basically expansion of (Od) /M d ~ 
(Aqcd/AO^ and is violated in small M 2 region. Thus we must take intermediate region 
where the integral is well-saturated by ground state contribution and the OPE is well- 
converged. These are quantified as: pole contribution > 50% of the total [5j |6j (for 
upperbound M^ ax ) and highest dimension terms < 10% of whole OPE (for lowerbound 
M^ in ). Practically, without these M 2 constraints, we will be stuck with the pseudo-peak 
artifacts leading to wrong conclusions, as discussed in [5J. 

Keeping these constraints in mind, let us consider how the first difficutly (i) is mani- 
fested in the QSR analyses. The multiquark correlators involve many quark lines, sharing 
only small fraction of the momentum with a scale relevant to the multiquark states. As a 
result, we have to incorpolate more non-perturbative diagrams in the OPE calculations 
than usual hadron analyses, to include relevant low energy correlations. Otherwise, we 
can not find the M 2 window due to the lack of sufficient low energy contributions. 

Here we shall add one remark about the low energy dominance. In the case of 
multiquark correlator analyses, it is frequently said that the pole-dominace is spoiled by 
large continuum (or high energy part of spectrum) contributions peculiar to multiquark 
correlators [BJ. This is a somewhat misleading expression. A real problem in the early 
studies is the lack of low energy correlation. In fact, in certain tetraquark correlators, 
inclusion of higher dimension terms enables to take sufficiently small M value, achieving 
the low energy dominance not much worse than usual meson sum rules, and even better 
than baryon sum rules [6j. Thus we conclude that a relevant problem is numerical 
ambiguities in values of higher dimension condensates. 

The second difficulty (ii) peculiar to multiquark cases arises as follows: even if we 
take into account the sufficient low energy correlations, they might be contributions from 
scattering backgrounds which are ground states in most of multiquark correlators. This 
prevents us from investigating the multiquark states of our interests. 

Both difficulties seem profound. However, as we will show, leading N c properties of 
multiquark states are tractable and, in principle, the corrections to them can be treated 
step by step with identifying origins of contaminations. 
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2 Factorization of condensates and separation of background based on 1/N C 



In this section, we discuss how useful 1 /N c expansion is in the application of the QSR. The 
first merit is presented in the evaluation of the higher dimension condensates. It is well- 
known that in large N c they can be factorized into the products of known condensates, 
(qq), (G 2 ), and (qg s aGq). Thus we can extend the OPE calculations to sufficiently higher 
dimension terms needed to include relevant low energy contributions. Here we must make 
one assumption on the N c scaling of these low dimension condensates. Since (O) =((qq), 
(a s G 2 ), {qg s crGq}) are known to be 0(N C ), we simply assume (0)|jv c = (0)N c /3. 

Next, we discuss how to separate the scattering background from the multiquark 
correlators. The point is that N c counting of the quark-gluon graphs is directly related 
to the qualitative picture of hadronic states. We will show that with taking appropriate 
interpolating fields, the multiquark state and background can be assigned in the different 
order of N c , and their mixing occurs only in the higher order of 1/N C than those studied 
in this work. Then we can concentrate on the quark-gluon graphs which are saturated 
by the isolated poles which reflect leading N c properties of multiquark states. 

To make discussions concrete, we consider the a meson with / = J — as a can- 
didate of the tetraquark (4q) state. The 4q operators with the a quantum number are 
given (assuming the ideal mixing for the a meson) by Jmm(x) — 53f=i ^m( x )Jm( x ) 
as products of meson operators = qrpTMq, where Dirac matrix Tm = (1)7/*) an d 
Tp (F = 1, 2, 3) are the Pauli matrices acting 0119= [u, d) T . 

Let us start the 1/N C classifications of quark-gluon graphs in the 2-point correlators, 
graphs (a)-(c) in Fig[TJ where planar-gluon lines are not explicitly drawn. The leading N c 
diagrams start from Fig[T]a), which include only 2 planar loops and thus reflect irrelevant 
free 2-meson scattering states. This indicates that the studies of 4q components require 
systematic steps beyond the leading N c . Thus we must proceed to the next leading order 
of 1/N C , 0(N C ) diagrams which could include the multiquark states. 
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Figure 1: 0{N 2 ) and 0(N C ) quark-gluon diagrams for 2 and 3 point correlators. 
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For systematic classifications, it is useful to first investigate the overlap strength of 
the operator J mm with hadronic states by 1/N C , then to use them in the classification 
of subleading diagrams of 2-point correlators. For this purpose, we consider the 3-point 
correlator among J mm and two separated Jm' (FigHJ d-f). Simple calculations indi- 
cate that the overlap strength of the 4q field with 2 meson states is (0| Jmm\M'M') = 
0(N c )S M m> +0(1) + .... since the leading order diagrams are 0(iV 2 ) (0(N C )) for M = M' 
(M ^ M') and the overlap strength of J' M with 2q meson state \M') is 0(iV c 1/2 ) [2]. 

Now we can classify the hadronic states in 2-point correlators (JmmJm'M 1 ) based on 
1/N C (See, Fig[l] a-c): (a) If M — M', O(N^) quark-gluon graphs include only the free 
2M scattering states in the region E > 2m m- If M ^ M', the contributions from these 
quark-gluon diagrams vanish, indicating the absence of free 2 meson scattering states, 
(b) 0(N C ) graphs include the 2M or 2M' scattering and could include 4q poles. 

Note that if M, M 1 arc different from pseudoscalar or axial vector, 2tt scattering 
intermediate states are not included up to 0(N C ) diagrams. Interactions are needed to 
transfer initial states into the 2tt intermediate states, but such interactions are suppressed 
by l/N c . Then the resonance peaks (if exist) below 2mM are isolated and have no width 
since the decay channels are absent at this order. Therefore, now we can reduce the a 
spectrum in the 4q correlator into peak(s) plus continuum if we retain only diagrams up 
to 0(N C ) for the appropriate currents. We will study these cases. 

This separate investigation of the O(N^) and 0{N C ) parts enables to perform the 
step by step analyses for each hadronic state. We relate the OPE, term by term of 1/N C , 
to the integral of the hadronic spectral function through the dispersion relation: 



1 ImITl„ (s) 

n^ e (-0 2 )= / ds- ^ (» = 2,1). (2) 



The ground state in the reduced 0(N C ) spectra Imlljv (s) can be described as the sharp 
peak due to the absence of the decay channels. Applying the usual quark-hadron duality 
to the higher excited states, 7rImII^ c (s) = A 2 <5(s — m\) + 9{s — st^)7rlmn^ e (s), and 
taking the moment of Borel transform of Eq. {2J , we can express the effective mass as 

mh(M ; Sth) = J^dse^i m u^( s) ■ (3) 

s t h can be uniquely fixed to satisfy the least sensitivity [4] of the expression ([3]) against 
the variation of M, since the physical peak should not depend on the artificial expansion 
parameter M. This criterion is justified only when the peak is very narrow, and our 
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1/N C reduction of spectra is essential for its application to allow the QSR framework to 
determine all physical parameters (raj, A, Sth) hi a self-contained way. 



3 Borel analyses for the reduced hadron spectra 



In this section, we show the results of Borel analyses. We concentrate on the 0(N C ) results 
of the off diagonal correlator (JvvJss)i whose leading order is 0(N C ) thus without the 
factorization violations at the 0(N C ) OPE. The OPE is carried out up to dimension 12. 
The numerical analyses are performed with the values for N c = 3 case, a s (lGeV) = 0.4, 
(a s G 2 /Tr) = (0.33 GeV) 4 , (qq) = -(0.25 ± 0.03 GeV) 3 , and to 2 , = (qg s aGq) / (qq) = 
(0.8 ± 0.1) GeV 2 , respectively. Most results shown below will use the central value. 

First we show in the left panel of Fig[2] the results of the large N c 2q correlators 
(expanded up to dimension 6) for the vector meson as a reference and the scalar meson 
as the 2q state in the a meson. The downarrow (upperarrow) indicates the values of 
M min ( M max( s th))- Following the E th (= ^fsE) fixing criterion, we fix E th to 1.0 (1.4) 
GeV for vector (scalar) mesons, and determine the mass as 0.65 (1.10) GeV. 

Now we turn to the 0(N C ) part of the 4q correlator, (JvvJss)> to investigate the 
possiblity of 4q states. Shown in the middle panel of Figf5] are the effective masses for 
E t h=1.0, 1.2, and 1.4 GeV. We select the £^ = 1.2 GeV case and evaluate the mass as 
~0.90 GeV, which is obviously lower than that of the 2q scalar meson case, ~1.10 GeV 
in large N c limit, and thus is considered as the mass of 4q state. Although these absolute 
values depend on the details of condensate values, the inequality m p < m^ q < ms is 
insensitive to them. In the right panel of Fig[2j we show the (qq) and to 2 , dependence of 
2q vector, scalar meson masses, and of the 4q mass deduced from (Jyv j\s)- 



Vector (S=1, L=0), 0(Nc| only 
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Figure 2: The large iV c 2q correlator results for the scalar and vector mesons (left), 
the 0(N C ) part of the 4q correlators (middle), and their mass relation for the various 
condensate values (right). 
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We have also investigated both 0(N 2 ) and 0(N C ) parts of (JssJss) ((^vvJvv))- As 
for 0(N 2 ) part, low energy contribution is almost absent below two meson thresholds. 
More precisely speaking, acceptable effective mass plots can be obtained only if we take 
much larger E t h value than the two meson thresholds, giving large effective masses. This 
indicates that 0(N 2 ) part includes only free scattering states consistently with 1/N C 
arguments. As for 0(N C ) part, (JssJss) ({•^vv^vv)) S* ve similar values ~ 0.9 (0.8) 
GeV, as (JvvJss) case i despite of the possible factorization violation in the formers. 

We have also investigated 4q correlators with 1 = 2. They do not allow us to find 
any stable effective mass plots or any reasonably wide M 2 window. This would indicate 
the absence of 4q pole in I = 2 channel, consistenly with experimental observations. 

In conclusion, it is found that the 1 /N c expansion is very useful to analyze the leading 
N c properties of multiquark states with separating scattering backgrounds. The QSR 
results show consistencies with expectations from 1/N C arguments, and further provide 
the prediction on 4q components in the a meson without 7T7t correlations. The leading N c 
analyses for multiquark correlators may be also useful to distinguish two meson molecules 
and tetraquarks, which is discussed for the recently found charmed mesons (X,Y,Z). Since 
initial states are controled by interpolating fields, and 1/N C countings give observations 
on intermediate states, the combination of them enables to separate the molecular type 
contributions in the spectral functions. Further discussions will be made elsewhere. 
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